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Abstract 



I The truncated Israel-Stewart theory of irreversible thermodynamics is 

I used to describe the bulk viscous pressure and the anisotropic stress in a 

' class of spatially homogeneous viscous fluid cosmological models. The gov- 



erning system of differential equations is written in terms of dimensionless 
variables and a set of dimensionless equations of state is utilized to complete 



the system. The resulting dynamical system is then analyzed using standard 

^ \ geometric techniques. It is found that the presence of anisotropic stress plays 

H 

I a dominant role in the evolution of the anisotropic models. In particular, 

in the case of the Bianchi type I models it is found that anisotropic stress 
leads to models that violate the weak energy condition and to the creation 
of a periodic orbit in some instances. The stability of the isotropic singular 
points is analyzed in the case with zero heat conduction; it is found that there 
are ranges of parameter values such that there exists an attracting isotropic 
Friedmann-Robertson- Walker model. In the case of zero anisotropic stress but 
with non-zero heat conduction the stability of the singular points is found to 
be the same as in the corresponding case with zero heat conduction; hence the 
presence of heat conduction does not apparently affect the global dynamics 



of the model. 
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I. INTRODUCTION 



Recently spatially homogeneous and isotropic imperfect fluid cosmological models were 
investigated using techniques from dynamical systems theory In these imperfect fluid 
models the bulk viscous pressure satisfies the truncated Israel-Stewart theory of irreversible 
thermodynamics. However, the models studied in |I| are isotropic Friedmann-Robertson- 
Walker (FRW) models which do not allow processes such as shear viscous stress; conse- 
quently the next step in this research programme is to study anisotropic generalizations of 
the FRW models. It was argued in [|I| that the anisotropic models studied in that allow 
shear viscous stress and which satisfy the Eckart theory of irreversible thermodynamics are 
not satisfactory since viscous signals in the fluid could travel faster than the speed of light - 
therefore it was concluded that a causal theory of irreversible thermodynamics such as the 
truncated Israel-Stewart theory ought to be utilized. 

Assuming that the universe can be modelled as a simple fluid and omitting certain 
divergence terms, the truncated Israel- Stewart equations for the bulk viscous pressure, 11, 
the heat conduction vector, qa, and the anisotropic stress, TTab (or shear viscous stress), are 
given by 0: 

n = -CKa + /5on-«oC)' 

= -^iTh^\T-^T,b + Ub + Piqb - aoU-b - aivr'^.J, (1.1b) 

Tlab = -2r]{Ua;b + P2Kb " aiqa;b) , (1-lc) 

where Ua is the fluid four-velocity, hab = UaUf, + gab is the projection tensor and (Aab) = 
^h'^a^'l{Acd + Adc — "^hcdh^^ Aef) . The variable T represents the temperature, k represents 
the thermal conductivity, /3o, j3i and (32 are proportional to the relaxation times, ao is a 
coupling parameter between the heat conduction and the bulk viscous pressure, and ai is a 
coupling parameter between the anisotropic stress and the heat conduction. We shall refer 
to equations ( |1.1| ) as the truncated Israel- Stewart equations. Equations ( |1.1| ) reduce to the 



Eckart equations used in B-H when ao = «i = /?o = A = = 0. 



Belinskii et al. were the first to study cosmological models using the truncated Israel- 
Stewart theory of irreversible thermodynamics to model the bulk viscous pressure and the 
anisotropic stress. Using qualitative analysis, Bianchi type I models were investigated sub- 
ject to equations ( p. . 1| ) in which equations of state of the form 

C = CoP'", r/ = 77oP", /9o = p-\ and = p-\ (1.2) 

were assumed , where m and n are constants and Co and rjo are parameters. The isotropiz- 
ing effect found in the Eckart models no longer necessarily occurred in the truncated Israel- 
Stewart models. It was also found that the cosmological singularity still exists but is of a 
new type, namely one with an accumulated "visco-elastic" energy 0. 

Recently, Romano and Pavon p,^] have studied anisotropic cosmological models in a 
causal theory of irreversible thermodynamics, analyzing the stability of the isotropic singular 
points in the Bianchi type I and III models. They also assumed equations of state of the form 



2|). However, they concluded that any initial anisotropy dies away rapidly but the shear 
viscous stress need not vanish and hence neither the de Sitter models nor the Friedmann 
models are attractors. 

In this paper we shall analyze qualitatively a class of anisotropic cosmological models 
arising from the use of the truncated Israel- Stewart equations, thereby expanding the anal- 
ysis in m to anisotropic models and extending the analysis in to causal theories. We 
will analyze both the Bianchi type V and the Bianchi type I models, which are simple gen- 
eralizations of the open and fiat FRW models. The system of ordinary differential equations 
governing the models is a dynamical system. We will find the singular points and determine 
their stability. In previous work dimensionless variables and a set of dimensionless 

equations of state were employed to analyze various spatially homogeneous imperfect fluid 
cosmological models. We shall utilize the same dimensionless variables and dimensionless 
equations of state here. One reason for using dimensionless equations of state is that the 
equilibrium points of the system of differential equations describing the spatially homoge- 



neous models will represent self-similar cosmological models [|l3]. In addition, it could be 



argued that the use of dimensionless equations of state is natural at least in some physical 
situations of interest (see e.g. fill)- (See also Coley |T^, and Coley and van den Hoogen 
for further motivation for using these dimensionless equations of state.) 

In section ^ we define the models and establish the resulting dynamical system. In 
section ^ we investigate the qualitative behaviour of the system for different values of the 
physical parameters. In particular, in section [III A| we analyze the system when there is no 
heat conduction and in section |III B| we analyze the system when there is heat conduction 



and bulk viscosity but with no anisotropic stress. In section |3 we discuss our results and 
make a few concluding remarks. For simplicity we have chosen units in which 8ttG = c = 1. 



II. ANISOTROPIC MODELS 

The Bianchi type V models are anisotropic generalizations of the negative curvature 
Friedmann- Robertson- Walker (FRW) cosmological models. The diagonal form of the Bianchi 
type V metric is given by: 

ds^ = -de + a{tfdx^ + hitfe^'^dy^ + c{tfe^''dz'' . (2.1) 

We assume that the fluid is moving orthogonal to the homogeneous spatial hyp ersurf aces; 
that is, the fluid 4- velocity, m", is equal to the unit normal of the spatial hypersurfaces. The 
energy-momentum tensor can be decomposed with respect to u"" according to [|1^]: 

Tab = {p+ P)UaUb + P9ab + QaUb + UaQb + T^ab, (2.2) 

where p = p + TI and p is the energy density, p is the thermodynamic pressure, 11 is the 
bulk viscous pressure, Hab is the anisotropic stress, and is the heat conduction vector as 
measured by an observer moving with the fluid. 

The Einstein field equations and the energy conservation equations in terms of the 
expansion(6') and shear (a) are: (see Coley and van den Hoogen 0): 

e = -2a^ -\d'' -]-{p + ?>p + m), (2.3a) 
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p = -e{p + p + n) - - 1 ( (Ti (2ni - u^) 

&i = -6ai + Hi, 
0-2 = -9a2 + 112, 



a2(2n2 - Hi; 



(2.3b) 
(2.3c) 
(2.3d) 



^2 = 3^2^3^^ 9^ (2.4a) 



a 



qi = -(Ti - CT2, (2.4b) 

where = |(cri + 0-2)^ — o-icr2- We have used the property that both a"^ = and n"-^ = 
to define new shear variables ai = a^i — (t\ ^"^^ ^2 = — cr^s and new anisotropic stress 
variables Hi = vr^j^ — tt\ ^"^^ ^2 = t^^i — ^r'^g to simplify the system. 



The evolution equations for 11, Hi, and 112 in this particular model, derived from ( p..lD , 
and using ( p^) , are then given by: 

e+'^ao{(Ti + a2){9^ -3cr^ -3p)y (2.5a) 
i-^«i(ai + a2)(^'-3a2-3p)), (2.5b) 
2 - (ai + (T2) {9' - 3a' - 3p)) . (2.5c) 



n = - 


n 

PoC 




ni = - 


Hi 


1 


2vP2 


~J2 


112 = - 


^2 


1 


2vP2 





Note that the heat conduction qi is completely determined by the shear via equation (|2.4bD ; 



thus equation ( |l.lb|) is not needed for the determination of the asymptotic behaviour of the 
models. 



Now the system of equations ( |2.3| ) and ( |2.5| ), define a dynamical system for the quantities 



(0, p, (Ti, (72, n. Hi, 112) = X of the form X = F(X). This system of equations is invariant 
under the mapping (see Coley and van den Hoogen P,p!0|) 



9 \9, (Ti AfTi, 0-2 Aa2, p A^p, 

p^AV, n^A^n, Hi-^A^ni, n2-^A2n2, 

C^AC, v^Xr], /5o^A-2/3o, (32^X-'(32 

ao A~^ao5 — ^ A~^ai, t — > X~H, 
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(2.6) 



and this invariance implies that there exists a symmetry in the dynamical system and 
hence a change of variables such that one of the equations can be made to decouple from 
the system. 

We define new dimensionless variables x, Si, E2, y, zi, Z2 and a new time variable Q as 
follows: 

Vsni v^n2 dQ 1 

-1-^, -2 = ^, and ^ = -/ (2.7) 

With this choice of variables, the Raychaudhuri equation, (p.3a| ), effectively decouples from 
the system. 

In order to complete the system of equations we need to specify equations of state for 
the quantities p, (, r], Pq, P2, Q^o and ai. In principle, equations of state can be derived 
from kinetic theory, but in practice one must specify phenomenological equations of state 
which may or may not have any physical foundations. Following Coley [|r^,|T^], we introduce 
dimensionless equations of state of the form: 

P_ _ I 

aia^'S 7^ = «2a;'■^ (2.8) 



^^-dio: , d2X , 

where po, Co, Voi ca and di {1 < i < 2) are positive constants, and i, m, n, and Pi {I < i < 2) 
are constant parameters {x is the dimensionless density parameter defined earlier). In the 
models under consideration 6 is strictly positive, thus equations (|2.8|) are well defined. 

We define new constants 61 = cti/Co, ^2 = 0,2/^0, Ci = aidi and C2 = a2'^2- The system of 
equations (|2.3| ), and (|2.5| ), written in the new dimensionless variables (|2.7|) and employing 
the above equations of state (|2.8| ), then become 

— =x{l- 2q) + 9pox' + y + Si(2;2i - Z2) + ^2(2^2 - ^1) 



^ 'Si + S2)(4-4a;-S2), (2.9a) 



4^3 

^ = Si(2-g)-12^i, (2.9b) 

^ = S2(2-g)- 122:2, (2.9c) 

^ = yibix''-"" - 2 - 2g) + 9aix''^ + CixP^+"i(Si + S2)(4 - 4x - S'), (2.9d) 

^ = 2i(262x'-^"" - 2 - 2g) + a2X^^Si - C2X^'2+''HSi + S2)(4 - 4a; - S^), (2.9e) 

^ = ^2(262a:^^-" - 2 - 2g) + 02X^^X3 - C2X^'^+'-HSi + S2)(4 - 4a; - S^), (2.9f) 

where = |(Si + S2)^ — SiS2- The quantity g is the generahzed dimensionless deceleration 
parameter given by 

<l = ^ = l{^ + y + ^Pox' + T?) , (2.10) 

where £ is the average length scale of the universe (i.e., 9 = 3|). 

Finally, from the Friedmann equation, ( p.4aj ), we obtain the inequality 



4-4x-S2 = 4^>0 (2.11) 

where equality implies that the model is of Bianchi type I. The interior of the parabola 
4 = + 4x in the phase space represents models of Bianchi type V, while the parabola 
itself represents models of Bianchi type I. There are other physical constraints that may be 
imposed, for example the energy conditions |T^, which may place bounds on the variables 



X, y, Si, S2, zi, and Z2. A full list of the energy conditions is given in Appendix 0. In the 
present work we shall always assume that x > 0, which states that the energy density in the 
rest frame of the matter is non-negative. This is a necessary condition for the fulfillment of 
the weak energy condition (WEC) |ll6| . 

The equilibrium points of the above system all represent self-similar cosmological models, 
except for those singular points that satisfy O/O"^ = — (g + l)/3 = 0. If g 7^ —1, the nature 
of the equations of state (W^ ) at the equilibrium points, is independent of the parameters /, 
m, n, ri, r2, pi, and p2, and is given by 
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/■ - - 

p oc p, Q (X p'2 , oc p2 



(2.12) 



Po (X p ^, P2 c<- P ^, oc p ^, and ai oc p ^. (2-13) 

Therefore natural choices for /, m, n, ri, r2, pi and p2 are respectively 1, 1/2, 1/2, 1, 1, — 1, 
— 1. We note that if there exists a singular point with q = —1, then it necessarily represents 
a de Sitter type solution which is not self-similar. 

The most commonly used equation of state for the pressure is the barotropic equation of 
state p = (7 — l)p, whence from ( |2.8|) po = |(7 — 1) and / = 1 (where 1 < 7 < 2 is necessary 



for local mechanical stability and for the speed of sound in the fluid to be no greater than 
the speed of light). In addition, / = 1 reflects the asymptotic behaviour of the equation of 
state for p. 

To further motivate the choice of the parameter ri, we consider the velocity of a viscous 
pulse in the fluid [0, 

where v = 1 corresponds to the speed of light. Using ( p. 7] ) and equations (|2.8D, we obtain 

V = {aix'^-^y/^. (2.15) 

Now, if ri = 1, then not only do we obtain the correct asymptotic behaviour of the equation 
of state for the quantity /?o but we also choose < ai < 1 since then the velocity of a viscous 
pulse will be less than the velocity of light for any value of the density parameter x. In this 
way the parameter ai has a physical interpretation as the square of the speed of a viscous 
pulse in the fluid. Therefore, in the remainder of this analysis we shall choose ri = r2 = 1. 
We shall also choose m = n = 1 and pi = p2 = —1 for simplicity. 

Using these particular values for m, n, ri, r2, pi, and p2, we can easily show that 
all singular points are self-similar except in the case 7 = 3(o, whence the singular point 
(x. Si, S2, y, -Zi, 2^2) = (1,0,0,-37,0,0) represents a de Sitter model. This is precisely the 
same as in the case in which the Eckart theory was employed 
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The full six- dimensional system is very difficult to analyze completely, so various phys- 
ically interesting subsystems are investigated. The case of zero heat conduction implies, 
via equations ( p.4b| ) and ( |2.7] ) that Si + S2 = 0. In addition, adding equations ( p.9b| ) and 



( p.9c| ) we deduce that zi + Z2 = 0, in which case the resulting system is four- dimensional 



see Section [111 Ap. The case of non-zero heat conduction with zero anisotropic stress is a 



three dimensional system and is discussed in Section [IIIB 



III. QUALITATIVE ANALYSIS 

The qualitative analysis of a system of autonomous ordinary differential equations X = 
F(X) begins with the determination of the singular points, that is, those points X such 
that F(X) = 0. To analyze the stability of each singular point we linearize the system in a 
neighborhood of each singular point. The signs of the eigenvalues of the derivative matrix 
determines the stability of the singular point, provided that the eigenvalues have non-zero 



real parts. (See Hirsch and Smale ||T8[, Sansone and Conti |]19|, Wiggins and Andronov 
et al. 1^ for useful reviews.) 

A. Zero Heat Conduction 

In the case of zero heat conduction, qi = 0, the field equations imply that Si + S2 = 
and Zi + Z2 = 0. Also, S^ = S^^^, hence (for simplicity), we shall drop the subscripts on S 
and Z] that is, S = Si = — S2 and z = zi = —Z2. The system of equations then becomes: 

x' = x(37 - 2 - 2q) + y + 6zE, (3. la) 

S' = S(2 - g) - 122, (3.1b) 

= y(6i - 2 - 2g) + 9aix, (3.1c) 

z' = z{2b2 -2-2q) + agxS, (3. Id) 

where 
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g = l((37-2)x + |/ + S2), (3.2) 

and the physical phase space is 

4 - 4a; - > and x>0. (3.3) 

There exists three obvious and physically motivated invariant sets in the phase space of 
the system. They are J^7^W := {{x,E,y, z)\J: = z = 0}, BI := {(x, S, ;z)|4 - 4x - = 
0, and S 7^ 0}, and BV := BI'^HJ-'TZW^ (where subscript c denotes the complement) which 
represents the Bianchi type V models. The set J^TZW represents the spatially homogeneous 
and isotropic negative and fiat curvature FRW models and the set BX represents the Bianchi 
type I models. There are possibly eleven different singular points of the system ( |3.1| ). The 
singular points lying in the set J-'TZW are 

(0,0,0,0), (l,0,i/-,0), (l,0,i/+,0), (3.4) 

where 

y^ = ^-^^±yih-3^y + 36a,. (3.5) 

Also, if Si = then there is a non-isolated line of singular points that passes through the 
points (0, 0, 0, 0) and (1, 0, y^, 0), where Bi is, 

Si = (37-2)(2-6i) + 9ai. (3.6) 

These points represent open (x = 0) and flat [x = 1) FRW models. There are possibly six 
singular points lying in the BI invariant set. They are 

(0,-2,0,0), (0,+2,0,0) (3.7) 

which represent Kasner models, and 

(x+, y+, +Z+), (x+, y+, -z+) 

(x",+S",?/",+z"), {x~,-f:~,y~,-z-) (3.8) 
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where 



and is given by 



where 



002 

t = f4a2 + (2 - 7)(&2 - 1 - t: 



2a2 



it -2), 
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(Si + 2S2) ± ^{Si - 2^2)2 + 96aia2 
4(2-7) ' 



(3.9) 



5i = (2-7)(6i-2) + 3ai, 
^2 = (2-7)(&2-l)+4a2. 

The remaining two singular points he in the BV invariant set. They are 



:i -60 



3a2 '^^3a2(2-6i)' 02(2-61) '^6^3a2(2-6i 
--62) 



(1-62)^1 3ai(l-62) 1 

' ''1 



(1 - h)B, 



-(1-62)51 3ai(l-62) 1 



1-62)^1 



(3.10) 
(3.11) 



(3.12) 



3a2 ' ^ 3a2(2- 61) ' 02(2-61) ' 6^18a2(2-6i) 

where Bi is given by equation 

The stabihty of each of these singular points is very difficult to determine in general. 
However, one question that can be asked is whether these anisotropic models generally 
isotropize; that is, "Does there exist a stable (t-time) equilibrium point in the set JFT^-W? 
(see subsection |III A 1|) . We shall also analyze the model when there is zero anisotropic 
stress {z = 0) in order to determine the effects that bulk viscous pressure may have on the 
models (see subsection [III A 21) . We shall also analyze the effect of anisotropic stress in an 
anisotropic model with zero bulk viscous pressure {y = 0). (see subsection [III A3| ). 
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1. Stability of Isotropic Singular Points 



In this subsection we are going to resolve the stabihty of the isotropic singular points, 
that is, those singular points lying in the set J-'TZW. We want to determine if there exists a 
stable (t-time) singular point in the future. In f2-time this translates to showing that there 
exists a source in the set J^TZW. 

The singular point (0, 0, 0, 0) represents the Milne model, and has eigenvalues 

2, 2(62 - 1), I {(37 + &i - 4) + ^^(37 + &i - 4)2 + 45i} , 

1 {(37 + 61 - 4) - ^(37 + &i-4)2 + 45i} , (3.13) 

where Bi is given by equation (|3.6|) . The bifurcation values are 62 = 1 and Bi = 0. If 
Bi = then there exists a non-isolated line of singular points. The stability of this point is 
summarized in Table |. 

The singular point (1, 0, y~, 0) represents a flat viscous fluid FRW model and has eigen- 
values 

~ {{bi + 37 - 4) - ^(6, + 37-4)2 + 4^1 1 , ^{b, + 3^ - 4)^ + 4B,, 

i {53 + v^^T^^} , ^{^3-^53'- 8^2}, (3.14) 

where 

B2 = (262 - 37 - y-)i6 - 37 - y-) + 24a2, (3.15) 
^3 = 462 - 3(37 - 2) - 3?/". (3.16) 

The stability of this point is summarized in Table |I|. 

The singular point (1, 0, y^, 0) represents a flat viscous fluid FRW model and has eigen- 
values 

_1 + 37 - 4) + v'(6i + 37-4)2 + 45i| , - ^'(fei + 37 - 4)2 + 45i, 

i {55 + ^^5^-854}, \[b,-^B,'-8B,Y (3.17) 
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where 

B, = (262 - 37 - 2/+)(6 - 37 - y+) + 24a2, (3.18) 
55 = 462 - 3(37 - 2) - 3y+. (3.19) 

The stabihty of this point is summarized in Table 

From the stabihty analysis of these singular points we can conclude that there exists a 
range of parameter values such that one of the singular points in the set J^TZW is a source 
(sink in t-time). If i?i < and 62 > 1 then the point (0,0,0,0) is a source — this result 
is similar to the observation of Coley and van den Hoogen [Q] in the Eckart theory when 
m = n = 1 and 9Co - (37 - 2) < 0. If Bi > 0, B2 > 0, and ^3 > then the point (1, 0, y-, 0) 
is a source. However, if either of these two conditions are not satisfied then the anisotropic 
models will not tend to an isotropic FRW model to the future (t-time) . Romano and Pavon 
remarked that the anisotropy dies away quickly in the anisotropic models and hence 
the cosmological model isotropizes, nonetheless the cosmological model does not tend to an 
FRW or de Sitter model. The same result is true here for some range of parameter values. 
If 62 < 1 and i?i < then the models all isotropize but the anisotropic stress does not tend 
to zero and therefore the model does not asymptotically approach an FRW model. 



2. Zero Anisotropic Stress 

In order to observe the effects of bulk viscous pressure in the model we set Hi = 112 = 0, 
which implies that 1/(32 = 0. In the model under consideration here this amounts to setting 
z = and a2 = in system (|3.1|) . The resulting system is three-dimensional and has the 
form 

x' = x{3-f -2-2q)+y, (3.20a) 
S' = S(2-g), (3.20b) 
y' = y{hi - 2 - 2g) + 9aix, (3.20c) 
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where 

g=i((37-2)x + y + S2), (3.21) 

and the physical phase space is 

4 - 4a; - > 0, and x> 0. (3.22) 

In this case there exist four invariant sets of particular interest. Similar to the previous 
analysis we have the set J^IZW := {(a;, S,?/)|S = 0} and BT := {(x, S,?/)|4 — 4a; — = 
and S ^ 0}. The Bianchi type V invariant set can be subdivided into two disjoint sets, 

i3V+ := BT n J^nw n {{x,j:,y)\j: > 0} and BV~ := BT n J^nw n {{x,i:,y)\j: < 0}. 

Due to the symmetry in the equations (reflection through the S = plane) the qualitative 
behaviour in the set S < is equivalent to that in the set S > 0; henceforth ( and without 
loss of generality), we shall only concern ourselves with the part of the phase space with 
S > 0. 

The singular points of the system (|3.20|) are 



(0,2,0), (0,0,0), (1,0,1/-), (1,0,1/+), (3.23) 



where ?/+ is given by equation 

There is only one singular point in the invariant set S > 0. The singular point (0, 2, 0) 
is in the set BT and has eigenvalues 



&1 + 37-12 1 r- ^ 

-4, -^(6i -37)2 + 36ai 

1 
2 



61 + 37-12 1 



+ -7(6i-37)2 + 36ai. (3.24) 



This singular point is either a saddle or a sink depending on the value of the parameter 

Eg = (2-7)(6i -6) + 3ai. (3.25) 

If Bq > 0, then the point is a saddle with a 2-dimensional stable manifold. If -Be < 0, 
then the point is a sink, and if = 0, then the point is degenerate (discussed later). The 
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solution at this singular point is a Kasner model. The stability of this point is summarized 
in Table [V[ 

The singular point (0, 0, 0) has eigenvalues 

2, ^ I (6i + 37 - 4) + v'(6i + 37-4)2 + 451 1 , 

^ + 37 - 4) - + 37-4)2 + 4^1 1 . (3.26) 

This point is either a saddle or a source depending on the value of the parameter Bi. If 
-Bi > 0, then the point is a saddle with a 1-dimensional stable manifold. If i?i < 0, then 
the point is a source, and if i?i = the point becomes degenerate (discussed later). The 
stability of this point is summarized in Table |V. 
The singular point (1,0, y^) has eigenvalues 

-I [{bi + 37 - 4) - ^(6i + 37-4)2 + 4^1 1 , .^(6, + 37 - 4)2 + 45i, 

I (61 + 37 - 12) - \/(6i-37)2 + 36ai| . (3.27) 

This point is either a saddle or a source depending on the value of the parameter Bi. If 
-Bi < 0, then the point is a saddle point with a 1-dimensional stable manifold. If Si > 0, 
then the point is a source, and if Si = the point becomes degenerate (discussed later). 



The stability of this point is summarized in Table |V. 
The singular point {1,0, y~^) has eigenvalues 

[ih + 37 - 4) + y^(6i + 37 - 4)2 + 4Si| , -^(61 + 37 - 4)2 + 4Si, 
-\ [(bi + 37 - 12) + ^Jih - 37)2 + 36ai| . (3.28) 

This singular point is either a saddle or a sink depending on the parameter Bq. li Bq < 0, 
then the point is a saddle with a 2-dimensional stable manifold. If Bq > 0, then the point 
is a sink, and if -Bg = then the point is degenerate (discussed later). The stability of this 
point is summarized in Table [rV| . 

The bifurcations in this model occur at -Bi = and Bq = 0. If Si = then there exists 
a line of singular points passing through the points (0, 0, 0) and (1, 0, y^). This line can be 
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shown to have some saddle-hke properties. In particular, if Si = then the points {1,0, y~) 
and (0, 0, 0) experience a saddle-node bifurcation. If i?6 = then the curve y = 3(2 — 7)0;, 
= 4 — 4x which lies in the set BI is singular. This observation is analogous to the case 
7 = 2 in perfect fluid Bianchi type V models 0. In particular, if i?6 = then the points 
(1,0,?/"'") and (0,2,0) experience a saddle-node bifurcation. 



All information about the singular points is summarized in Table [I^. It is very easily 
seen that if -Be > then the solutions tend to an isotropic model both to the past and 
to the future (in t-time), while if -Be < 0, then solutions only tend to an isotropic model 
to the future. Provided that there are no periodic or closed orbits in the set S > 0, all 
models isotropize to the future {Q —00 or t ^ 00). Note the difference in the result 
here and the result in the previous subsection. If there is a 'non-zero' anisotropic stress 
then there is a range of parameter values such that models will not isotropize, and if there 
is 'zero' anisotropic stress then all models will isotropize to the future. Therefore we can 
conclude that in the truncated Israel- Stewart theory the anisotropic stress plays a dominant 
role in determining the future evolution of the anisotropic models. This result is contrary to 
the observations in Coley and van den Hoogen based upon the Eckart theory where the 
anisotropic stress played only a minor role and did not determine the the future evolution 
of the models. 

3. Zero Bulk Viscous Pressure 



As we have seen in subsection |III A 1| , anisotropic stress plays a dominant role in the 



evolution of the anisotropic cosmological models. To further analyze the effects of anisotropic 
stress on the evolution of an anisotropic model we shall set 11 = 0, which in our model implies 



that l/Po = 0. This translates into setting y = and ai = in equations ( |3.1| ). In order 
to illustrate the possible influence anisotropic stress may have on an anisotropic model we 
further restrict ourselves to the set BX := {(x, S, z)\4: — 4a; — = 0}. The resulting system 
is planar and lends itself easily to a complete qualitative analysis. 



17 



Consequently the system under consideration is 

E' = ^(2-7)E(4-E2)-12^, (3.29a) 
8 

z' = (262 - 37)^ - ^(2 - 7)^E2 + ^E(4 - E"). (3.29b) 
The singular points are (0,0), (+2,0), (-2,0), (E+,z+), and {^'^z') where 



(3.30) 

where 

(2-7)(62-3) + 4a2. (3.31) 
The point (0, 0) represents a fiat FRW model; the eigenvalues of this point are 



\}^Bs±^B^- 48[B, + ^(2 - 7)2] I , (3.32) 

where 

S8 = 4&2 - 3(37- 2). (3.33) 

If ^7 + 1(2 - 7)2 < 0, then the point (0, 0) is a saddle point. If S7 + |(2 - 7)^ > 0, then 
the stabihty of the point (0, 0) depends on the parameter Eg. If -Bg > the point (0, 0) 
is a source and if < the point (0, 0) is a sink. Bifurcations of this point occur when 
i?7 = — 1(2 — 7)2 and i^g = and are discussed later. 

The points (±2, 0) represent Kasner models. The eigenvalues are 

^ I (262 + 37 - 12) ± ^(262 + 37 - 12)2 + 24^7 1 . (3.34) 

If i?7 > 0, the points (±2, 0) are saddle points. If B-^ < 0, then the points (±2, 0) are sinks. 
The bifurcation that occurs when S7 = is discussed later. 

The points (S^, z"^) only exist when S7 + |(2 — 7)^ > 0. The eigenvalues are 

-I I 262 + 37 - 12 - , ^ R7I 
2\\ (2-7) V 



5 / 2 
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If -|(2 - < ^7 < 0, then the points z^) are saddle points. If Bj > 0, then the 
points (S^*", z"^) are sinks. The bifurcation vakies Bj = — 1(2 — 7)^ and Bj = are discussed 
later. The stability of all singular points is summarized in Table |V]. 

Knowing the singular points and their eigenvalues only reveals the local behaviour of the 
system ( |3.29| ). The determination of some of the global properties requires investigating the 
existence or non-existence of periodic orbits and analyzing points at infinity. 

The existence of periodic orbits is difficult to prove. However, with the aid of Dulac's cri- 
terion |2l[|, we are able to prove the non-existence of periodic orbits for a range of parameter 
values. Taking the divergence of the system ( p.29|) we can see that 

V-/ = ^58-f (2-7)Sl (3.36) 

Therefore, if -Bs < then there do not exist any periodic orbits. 

To analyze the points at infinity we first change to polar coordinates = S^ + z^ and 6 = 
tan~^(2;/S) and then we compactify the phase space through the following transformations: 

f=^, -9 = 9, f = (l-r)^ (3.37) 

in which case the evolution equations for f and 9 become: 

df (3 - - - -1 

_ = (1 - f)¥ |-(2 - 7) cos^ 9 + (262 - 37) sin^ 9 + {a2 - 12) cos^ sin^j 

cos^ 9(3 - - - -1 

— (1 — f)f —— -(2-7) cos^ ^ + 3(2 - 7) sin^ 9 + a2Cos^sin^>, (3.38a) 

I (262 - ^7 - 3) cos ^ sin ^ + 03 cos^ 9 + 12 sin^ ^| 



d9 ,^ .0 { 3 



-f'^^^^(2-7)sin^ + a2CoseL (3.38b) 



4 12 

The points at r = 00 are mapped to the unit circle f = 1. Hence the singular points at 
infinity are those points on the unit circle f = 1 where ^ = 0. The singular points are thus 

(1,^), (1,-^), (l,r), and (l,r + 7r), (3.39) 

where tan 6'* = — 2a2/3(2 — 7). In order to determine the stability of these singular points 
we need to study the values of ^ and ^ in a neighborhood of each of the singular points. 
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We find that the points (1, 9*) and (1, 9* + vr) are saddle-points while the points (1, ±vr/2) 
are sources. 

To obtain a complete picture of the qualitative behaviour of the model we must discuss 
the various bifurcations that occur. A bifurcation occurs at Bj = — 1(2 — 7)^, in which 



case the point (0, 0) undergoes a pitchfork bifurcation |^0[ to create the two new singular 
points and its stability is transferred to them. When Bj > — 1(2 — 7)^, the point 

(0, 0) experiences an Andronov-Hopf bifurcation at -Bg = |2^. Therefore, it can be shown 
that there exists a 5 > such that for every Bg G (0, 6) there exists a periodic orbit. In 
addition, the periodic orbit is an attractor. A third bifurcation occurs at i?7 = when the 
points and (±2,0) undergo a transcritical bifurcation in which they exchange 

stability. The stability of all of the singular points (finite and infinite) is given in Table 0. 

Let us now discuss the qualitative properties of this model. If Bj < — 1(2 — 7)^, then all 
trajectories evolve from the singular points at (±2, 0) representing Kasner models to points 
at infinity. These models are generally unsatisfactory since the WEC (which implies < 4 
for the Bianchi I models here), is broken eventually. However, there do exist two exceptional 
trajectories for which the WEC is satisfied always. These are the trajectories that lie on the 
unstable manifold of the singular point (0, 0) which describe models that have a Kasner-like 
behaviour in the past (t-time) and isotropize to the future toward the point (0, 0). A phase 
portrait of this model is given in Figure ^ 

If — I (2 — 7)^ < Bj < and Bs < then there exist two classes of generic behaviour. 
One class of trajectories evolve from the isotropic singular point (0, 0) and evolve to points 
at infinity. The second class of trajectories evolve from the singular points (±2,0), which 
represent Kasner models, and evolve to points at infinity. Both of these classes of trajectories 
describe models that fail to isotropize, and describe models that will eventually violate the 
WEC. If — 1(2 — 7)^ < i?7 < and Bs < 0, then there exists three sets of exceptional 
trajectories. One set is the stable manifolds of the points (S^,z^) which represent models 
that start at (S^, z^) and evolve to points at infinity and hence the WEC will be violated. 
There do exist trajectories describing models that satisfy the WEC for all time, namely the 
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unstable manifolds of the point (E^,z^). One set of these trajectories start at the singular 
point (±2,0) which represent the Kasner models and evolve to the point (S^,^;^). The 
other set of trajectories evolve from the isotropic singular point (0, 0) to the singular points 
(S^, z^). In this case there are no models which isotropize. A phase portrait of this model 
is given in Figure ^ 

If — 1(2 — 7)^ < i?7 < and > then there exist three classes of models. One class of 
trajectories evolves from the periodic orbit to the isotropic singular point (0, 0). This class of 
models is interesting in that the past singularity has an oscillatory nature, that is, both the 
dimensionless shear S and the dimensionless anisotropic stress z tend to a closed periodic 
orbit in the past (t-time). This class of trajectories also represent models that isotropize 
and represent models in which the WEC is satisfied always. The second class of trajectories 
are those which evolve from the periodic orbit to points at infinity. This class of trajectories 
represent models that will not satisfy the WEC at some point in the future. The third class 
of trajectories is the same as the second class of trajectories in the case — 1(2 — 7)^ < Bj < 0. 
Again there exist three sets of exceptional trajectories. The stable manifolds of the points 
(S^,z^) represent models that will eventually violate the WEC. The unstable manifolds 
of the points represent either models that start at the Kasner-like singular-point 

(±2,0) and evolve to the point (S^,2;^) or represent models that start from the periodic 
orbit and evolve to the point {J]'^,z'^). The phase portrait in this case is very similar to 
that of Figure 

If -B7 > and Bs < 0, then the behaviour of the trajectories is very similar to the 
behaviour of the trajectories in the case — 1(2 — 7)^ < Bj < and Bs < 0. The difference 
stems from the fact that the points (±2, 0) are now saddles and the points (S^, z^) are now 
sinks. The phase portrait in this case is very similar to that of Figure |^. 

If i?7 > and B^ > 0, then the behaviour of the trajectories is very similar to the 
behaviour of the trajectories in the case — 1(2 — 7)^ < Bj < and Bs > 0. The difference 
stems from the fact that the points (±2,0) are now saddles and {'E'^,z'^) are now sinks. A 
phase portrait of this model is given in Figure |^. 
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In conclusion, the general behaviour of these models is unsatisfactory in that the WEC 
is violated eventually, except in the case -Bg > and — 1(2 — 7)^ < B7 where there exists a 
set of models (of non-zero measure) that will always satisfy the WEC. These are the models 
represented by the trajectories which start at the periodic orbit and isotropize to the point 
(0,0) to the future (t-time). There also exist models which satisfy the WEC always, but 
these are the models represented by the unstable manifolds of the saddle-points. 

Clearly the anisotropic stress in the truncated-Israel- Stewart theory plays a very domi- 
nant role in the evolution of the anisotropic models. This is in contrast to what was found 
in using the Eckart theory where it was found that anisotropic stress played a very minor 
role in determining the asymptotic behaviour. However, if -Bg < 0, then all models are 
generically unsatisfactory in that the WEC will be violated. If -Bg > 0, then there does exist 
a set of satisfactory models where the WEC will always be satisfied. It is also interesting to 
note briefly the existence of a periodic orbit, this type of behaviour is not seen in the Eckart 
models. With the existence of this periodic orbit, the past attractor which this periodic orbit 
represents, has a oscillatory character to it, in that the dimensionless shear (and therefore 
the dimensionless density) and the dimensionless anisotropic stress will have an oscillatory 
nature. 



B. Non-Zero Heat Conduction 

In this section we will study the effects of heat conduction on the models. For simplicity 
we will assume that the anisotropic stress is zero. Although Hab = 0, the bulk viscosity 
is still present, and hence, in a sense, we are investigating the effect heat conduction will 
have on the viscous models with bulk viscosity. Under these assumptions, the system ( p.9|) 
reduces to a one-parameter family (that is, in addition to the parameters arising in the 
equations of state) of three-dimensional systems 0. We label this new parameter C, which 
is a function of the integration constant that appears when equations ( p.9b| ) and ( |2.9c| ) are 



integrated. The parameter C = {k + l)/vP — ^H-T is bounded between — 1 < C < 2 and 
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is analogous to the parameter used in . The value C = corresponds to the case in which 
there is zero heat conduction, and C = 2 corresponds to the case in which the model is 
locally rotationally symmetric (LRS). We define a new shear variable 

whence the system becomes 

x' = x(37 - 2 - 2g) + ?/ - ^^(4 - 4x - S^), (3.41a) 
S' = S(2-g), (3.41b) 
y' = y{hi -2-2q) + 9aix + V3CciS(4 -Ax- T?), (3.41c) 



g=i((37-2)x + y + S2), (3.42) 



where 



2 

and the physical phase space is 

4 - 4a; - > and a; > 0. (3.43) 

There exists three physically interesting invariant sets in the phase space of the system, 
namely, mW := {(x,S,y)|S = 0}, BX := {{x,T.,y)\A - Ax - = 0, and S ^ 0}, and 
BV := BX'^ n J^TZW^ (where subscript c denotes the complement) which represents Bianchi 
type V models. As before, the set J-'TZW represents the spatially homogeneous and isotropic 
negative and fiat curvature FRW models and the set BX represents the Bianchi type I models. 

There are six different singular points of the system. The singular points lying in the set 
J^nW are: 

(0,0,0), (1,0,2/-), (1,0,1/+), (3.44) 



where is given by equation ( |3.5| ). Also, if i?i = then there is a non-isolated singular line 
that passes through the points (0,0,0) and (1,0,?/^), where Bi is given by equation (^.6[). 
These points represent open {x = 0) and flat (x = 1) FRW models. The eigenvalues of the 
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linearization in a neighborhood of each of the isotropic singular points are the same as in the 
case with zero heat conduction (see equations ( p.26| ), ( |3.27| ), and ( p. 281 ) for the eigenvalues 



of (0, 0, 0), (1, 0,y ), and (1, 0, y^), respectively, and the appropriate parts of Table |rV|). 
The singular points located in the set BI are 

(0,-2,0), and (0,+2,0). (3.45) 

The eigenvalues of the linearization about the point (0, —2, 0) are similar to those in the case 



with zero heat conduction [see equation (|3.24| )]; indeed, only the first eigenvalue is different. 



namely, instead of Ai = —4, we now have Ai = —4 — 2C, which is very easily seen to be 
negative definite. Therefore the stability of the point (0, —2, 0) is the same as in the previous 
case with zero heat conduction. Similarly, for the eigenvalues of the linearization about the 
point (0, +2,0), only the first eigenvalue is different, namely, instead of Ai = —4 we now 
have Ai = —4 + 2C, which is negative definite for C ^ 2. Therefore, if C 7^ 2, the stability 
of the point (0, +2, 0) is the same as in the case with zero heat conduction. The case C = 2 
is discussed below. 

The sixth singular point is {x, S,|/), where 

4(C2-4)(6i -6 + 4v^Ci) 



X 



C2[16V3ci + (37 - 2) (61 - 6) - 9ai] ' 
4 



12(C^-4)[4v/3ci(2-7)-3ai] 
C2[16v^ci + (37 - 2){h - 6) - 9ai] ' ^ ' ' 

However this last singular point lies outside the region of phase space defined by equations 



( ra ) for C^2. 

If C = 2, then a 'transcritical' bifurcation occurs, the points (x,S,?/) and (0, +2,0) 
coalesce and become a single point. The stability of this point cannot be determined via 
linearization. If (37 — 2) (6 — bi) + 9ai — 16-\/3ci = 0, then there is a line of singular 
points y + (37 — 2)x = 0, S = 2. The stability of the singular point is very difficult 
to determine analytically (even with the use of center manifold theory |2^). However, 
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numerical experiments in addition to some analysis show that the singular point has some 
of the same behaviour as in the case with C ^ 2 (e.g., if i?6 > the point is a saddle and if 
i?6 < then the point has both saddle-like and sink-like behaviour). 

The stability of the singular points (0, 0, 0), (1, 0, y'), (1, 0, ?/+), (0, -2, 0), and (0, +2, 0) 
are the same as in the case with zero heat conduction. The heat conduction does not 
determine the stability of the singular points that lie in the physical phase space ( p.43| ) but 
does play a role in determining their eigendirections. This is similar to the situation in which 
the bulk viscous pressure is absent whence the model reduces to one that was analyzed in 
Coley and van den Hoogen 0; there the addition of heat conduction did not change the 
stability of the singular points but did allow the models to violate the WEC by rotating the 
principal eigendirections. 



IV. CONCLUSION 



This work improves over previous work on viscous cosmology using the non-causal 



first-order thermodynamics of Eckart ||22| in that a causal theory of irreversible thermo- 
dynamics has been utilized. Also, this work enhances the analysis of anisotropic viscous 
cosmologies in 10-^] because more than just the isotropic singular points have been ana- 
lyzed. The present work also generalizes the analysis of causal viscous FRW models in 

i- 

Again we have seen that the singular points of the dynamical system describing the 



evolution of an anisotropic viscous fluid cosmological model are, in general, self-similar |]I0 
In the case in which m = n = 1, rj = 1 and Pi = —1 {i = 1,2) all singular points are 
self-similar except in the case in which 7 = 3(o when there exists a singular point that 
represents a de Sitter model which is not self-similar. 

We have found that in the case of zero heat conduction the anisotropic models need 
not isotropize (that is, there exists a range of parameter values and initial conditions such 
that the models will not isotropize). The parameter 62, which is the parameter related to 
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the relaxation time of the anisotropic stress, plays a major role in determining the stability 
of the isotropic models. In the special case of zero anisotropic stress we have shown that 
all models isotropize. The addition of anisotropic stress on an anisotropic Bianchi type 
I model reveals some of the effects that anisotropic stress has on an anisotropic model. 
For instance, anisotropic stress generically causes models to increase their anisotropy and 
eventually violate the weak energy condition. Anisotropic stress in some instances causes 
the creation of a periodic orbit. This periodic orbit represents a past attractor in which 
the dimensionless quantities p/O"^ and a/9 are approximately periodic. It is interesting to 
note that it is only when this periodic orbit is present that there exist trajectories (in the 
interior) which represent models that will isotropize and satisfy the weak energy condition. 

The models with heat conduction analyzed here had no anisotropic stress but did have 
bulk viscosity. Consequently, we have investigated whether any qualitative changes arise 
from the inclusion of heat conduction. From our analysis the addition of heat conduction 
in the model did not change the stability of the singular points, an hence the asymptotic 
states of the models. However, the inclusion of heat conduction did affect the dynamics in 
a neighborhood of each of the singular points since the eigendirections changed. 

In this work we have employed the truncated Israel- Stewart theory which is a causal and 
stable second order relativistic theory of irreversible thermodynamics. It is possible that the 
truncated theory is applicable in the very early universe. However, it is known that such 
a truncated theory could result in some pathological behaviour, (e.g., in the temperature 
) . Hence this work should be considered as a first step in the analysis of the full theory of 



Israel-Stewart-Hiscock p4| , p!7| ,|6|,p5| . The present paper and [|T| provide a firm foundation for 
the analysis of viscous cosmological models using the full theory which we hope to present 
in the future. 
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APPENDIX A: ENERGY CONDITIONS 



For an imperfect fluid the energy conditions can be formulated with respect to the 



eigenvalues of the energy momentum tensor [0. The weak energy condition (WEG) states 
that TabW"'W^ > for any timelike vector W"' |15|. In the models under investigation the 



WEG, written in dimensionless variables, becomes 

3(2 --f)x-y- 2V3{zi + z^) + 9A > 0, 
37a; + ?/ + 2^3(^1 - 5^2) + 9A > 0, 
37X + ?/ + 2^3(^2 - 5^i) + 9A > 0, 



(Al) 
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where 



A = 4Vi67 + 2?/ + eVsizi + Z2W - 3(Si + S2)2(4 - 4x - S2). (A2) 

io 



We note that the eigenvalues of the energy momentum tensor must be real |]T6[ and therefore 
the quantity under the square root sign in ( |7^ ) must be positive. 

The dominant energy condition (DEC) states that for every timelike W"", TabW°'W^ > 
and TIW'^ is non-spacelike |]15|. In the models under investigation the DEC becomes 



< 3(2 - 7)3; - y - 2^3(^1 + Z2) 

< 37X + y + 273(^1 - 5^2) + 9A < 6(2 - 7)x - 2y - 4^3(^1 + Z2) + 18A, 

< 37X + ?/ + 2^3(^2 - 5zi) + 9A < 6(2 - 7)^ - 2y - aV3{zi + Z2) + ISA. (A3) 



The strong energy condition (SEC) states that TabWW' - \T^^W^Wb > for any 



timelike vector W"" |15]. In the models under investigation the SEC becomes 



WEC and 6(7 - l)x + 2y ~ 2^3(^1 + Z2) + 9A > 0. 



(A4) 
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TABLES 

TABLE I. Stability of the singular point (0,0,0,0) where dim(VF'') is the dimension of the 



stable manifold with respect to Q-time. 

sgn(Bi) sgn(b2 — 1) type dim(VF*) 

+ + Saddle 1 

+ - Saddle 2 

— + Source 

Saddle 1 



TABLE IL Stability of the singular point (1, 0,y ,0) where dim(VF*) is the dimension of the 



stable manifold with respect to fi-time. 

sgn(Bi) sgn(B2) sgn(B3) type dim(l^*) 

+ - Saddle 1 

+ + + Source 

Saddle 2 

- + + Saddle 1 

- + - Saddle 3 



TABLE in. Stability of the singular point (1, 0, y'*', 0) where dim(l^'') is the dimension of the 
stable manifold with respect to Jl-time. 

sgn(B4) sgn(B5) type dim(Ty^) 

Saddle 3 

+ + Saddle 2 

+ - Sink 4 
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TABLE IV. Stability of the 


singular points with respect to J7-time (both with and witout heat 


conduction). 










sgn(B6) sgn(Bi) 


(0,0,0) 


(l,0,y- 


) 


(1,0,2/+) (0,2,0) 


+ 


source 


saddle 




sink saddle 


+ + 


saddle 


source 




sink saddle 




source 


saddle 




saddle sink 


+ 


saddle 


source 




saddle sink 


TABLE V. Stability of the 


singular points 


with respect to f2-time 


, at both finite and infinite 


values, for the Bianchi type I anisotropic model with anisotropic stress and zero viscous pressure. 


B7 sgn(B8) 


(0,0) (±2,0) 


(S±,^±) (0,±oo)^ 


(±00, ^poo)*^ periodic orbit 


57 < ^(2 -7)2 


saddle sink 




source 


saddle 


^(2-7)2<57<0 - 


sink sink 


saddle 


source 


saddle 


^(2 - 7)2 < < + 


source sink 


saddle 


source 


saddle sink 


^7 > 


sink saddle 


sink 


source 


saddle 


Br>0 + 


source saddle 


sink 


source 


saddle sink 


'^These are the points at infinity corresponding 


to (f = 1,^ = 


= ±7t/2). 





'^These are the points at infinity corresponding to (r = 1, ^ = 9*) and {r = 1,9 = 9* + tt). 
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FIGURES 

z 




FIG. 1. The phase portrait describes the quahtative behavior of the anisotropic Bianchi type 
I models with anisotropic stress and zero bulk viscous pressure in the case < — 1(2 — 7)^. The 
arrows in the figure denote increasing O-time (O — > 00) or decreasing t-time (t ^ O"*"). 
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z 




FIG. 2. The phase portrait describes the quahtative behavior of the anisotropic Bianchi type 
I models with anisotropic stress and zero bulk viscous pressure in the case — 1(2 — 7)^ < By < 
and Bg < 0. The arrows in the figure denote increasing Jl-time (f2 00) or decreasing t-time 
(i^0+). 
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FIG. 3. The phase portrait describes the quahtative behavior of the anisotropic Bianchi type I 
model with anisotropic stress and zero bulk viscous pressure in the case B7 > and Bs > 0. The 
closed elliptical orbit close to the center represents the periodic orbit. The arrows in the figure 
denote increasing 0-time (O 00) or decreasing i-time {t — > 0"'"). 
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